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Using pulsed optical excitation and read-out along with single phonon counting techniques, we 
measure the transient back-action, heating, and damping dynamics of a nanoscale silicon optome¬ 
chanical crystal cavity mounted in a dilution refrigerator at a base temperature of Tf ss 11 mK. 

In addition to observing a slow (~ 740 ns) turn-on time for the optical-absorption-induced hot 
phonon bath, we measure for the 5.6 GHz ‘breathing’ acoustic mode of the cavity an initial phonon 
occupancy as low as (n) = 0.021 ± 0.007 (Tf, ~ 70 mK) and an intrinsic mechanical decay rate of 
70 = 328 ± 14 Hz (Q m « 1.7 x 10'). These measurements demonstrate the feasibility of using short 
pulsed measurements for a variety of quantum optomechanical applications despite the presence of 
steady-state optical heating. 

PACS numbers: 42.50.Wk, 42.65.—k, 62.25.—g 

The recent cooling of nanomechanical resonators to their motional quantum ground state m opens the possibility 
of utilizing engineered mechanical systems strongly coupled to optical or microwave fields for a variety of quantum 
metrology and information processing applications [4j| , amongst them the preparation of highly non-classical mechan¬ 
ical states and coherent frequency conversion between microwave and optical signals [SHU]- A particularly 
interesting device architecture for realizing large radiation pressure coupling between light and mechanics is the thin- 
film optomechanical crystal (OMC) [T3J[2], in which optical and acoustic waves can be guided and co-localized via 
patterning of the surface layer of a microchip. Based largely upon the OMC concept, new ideas for phononic quantum 
networks m and optomechanical metamaterials m have been proposed, in which arrays of cavity-optomechanical 
resonators are coupled together via optical or acoustic degrees of freedom, and in which laser light is used to para¬ 
metrically control the emergent network or material properties. 

For many of the above mentioned applications, operation at cryogenic temperatures is desired as it offers a route 
to obtaining the requisite low thermal occupancy and long mechanical coherence time. Recent measurements at 
millikelvin (mK) bath temperatures of an OMC resonator formed from single crystal silicon (131 117] , however, have 
shown substantial mechanical mode heating and mechanical damping due to weak sub-bandgap optical absorption [lSj . 
Although optical Q-factors in excess of 10 6 are realized in these highly optimized structures m, the large impact of 
even very weak optical absorption can be attributed to a combination of the relatively large energy per photon, and the 
sharp drop in thermal (phonon) conductance with temperature in the low temperature limit. Further complications 
arise from the seemingly contradictory requirements of isolating the mechanical resonator from its environment to 
obtain high mechanical Q-factor, and that of providing large thermal anchoring to a low temperature bath for cooling 
of the mechanical resonator. 

In this work we utilize pulsed optical excitation and single phonon counting |19j to study the transient dynamics 
of optical back-action, heating, and damping of the 5.6 GHz mechanical mode of a silicon optomechanical crystal 
resonator at mK bath temperatures. Phonon counting, realized by photon counting of the optically filtered motional 
sidebands of the reflected optical excitation pulse, yields simultaneously a high time resolution (~ 25 ns) and me¬ 
chanical mode occupancy sensitivity (< 10” 2 ). Measurement of both Stokes and anti-Stokes sidebands also yields an 
absolute calibration of the occupancy of the resonator mode in terms of mechanical vacuum noise I20fl22| . In addition 
to measuring initial phonon mode occupancies as low as (n) = 0.021 ± 0.007 and mechanical decay times as long as 
r = 475 ±21 /zs, we observe a slow (~ 740 ns) turn-on time for the optical-absorption-induced phonon bath that both 
heats and damps the mechanical resonator mode. Taken together, these measurements demonstrate the feasibility of 
using short pulsed measurements for quantum optical state engineering of the mechanics in silicon optomechanical 
crystals, despite the presence of large steady-state optical heating. 
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FIG. 1. a, Pulsed pump light at frequency wi (red arrows) is directed to an optomechanical crystal cavity (OMC) inside a 
dilution fridge via an optical circulator. The cavity reflection is then filtered at the cavity frequency uj c (black arrows) and 
directed to a single photon detector (SPD). b, SEM image of the silicon OMC cavity studied in this work, c, Finite-element- 
method simulation of the localized acoustic resonance at 5.6 GHz of the OMC cavity. Deformation of the beam structure is 
exaggerated to highlight the mechanical motion, with color indicating the regions of high (red) and low (blue) displacement 
magnitude. The inset shows the three processes affecting the mode occupancy: (i) detuned pump light at u>i (red) exchanges 
energy with the acoustic mode at rate 7om, generating scattered photons at u) c (black) in the process, (ii) pump light drives 
excitation of electronic defect states in the silicon device layer, subsequently exciting a hot phonon bath which heats the 
localized acoustic resonance at rate 7 P , and (iii) phonons escape the cavity volume via the acoustic shield at intrinsic decay rate 
70 , coupling the localized resonance to the fridge bath, d, Noise-equivalent phonon number unep (gray squares) and phonon 
occupancy (n) (red circles) for red-detuned (A = w m ) continuous-wave (CW) pumping versus intracavity photon number n c . 
The red dashed line indicates the expected hnep contribution from SPD dark counts. The vertical dashed line at n c ~ 45 
(4.5 x 10 -5 ) indicates the photon number during the on-state (off-state) of the pulse. Solid green and purple lines show fits to 
the CW heating model for fridge bath temperatures of T{ = 70 mK and 10 mK, respectively. 


The device studied in this work consists of a patterned nanobeam, formed from the top silicon device layer of a 
silicon-on-insulator wafer. The etched hole pattern in the silicon nanobeam forms an optomechanical crystal, in which 
photonic and phononic bandgaps at the ends of the beam co-localize optical and acoustic resonances with frequencies 
of u! c /2tt ss 196 THz (free space wavelength A c « 1530 nm) and w m /27r ss 5.6 GHz, respectively. A scanning electron 
micrograph (SEM) of the device is shown in Fig. [TJa, and a finite-element method simulation of the acoustic resonance 
is displayed in Fig. [TJ:. Coupling to the optical resonance is accomplished via an end-fire coupling scheme, using a 
lensed optical fiber inside of a dilution refrigerator to couple to the central waveguide shown in the SEM, as described 
in Ref. (18]. Surrounding the cavity is a 2D cross pattern [14] which possesses a complete acoustic bandgap in the 
5 — 6 GHz range, providing an additional acoustic shield for the mechanical resonator mode while allowing phonons 
above and below the phononic bandgap to carry heat from the nanobeam structure. 

The experimental set-up is shown schematically Fig. ID 1 - Electro-optic modulation of a laser probe beam generates 
optical pulses with frequency and duty-cycle controlled by a variable delay electrical pulse generator. The red arrows 
indicate a coherent pump at frequency uq which is red-detuned from the optical resonance frequency « c by A e 
u> c — uq = uj m . In this case absorption of a single phonon from the nanomechanical resonator results in upconversion 
of a pump photon to the anti-Stokes sideband at ui c , represented by the black arrows. The cavity reflection is filtered 
to reject the pump frequency and subsequently directed to a high-efficiency single-photon detector (SPD) and a time- 
correlated single photon counting system synced to the pulse generator. This measurement repeats each pulse period, 
building up a histogram with respect to photon arrival time relative to the sync pulse during a certain integration time. 
As the number of anti-Stokes photons is directly proportional to the average phonon occupancy of the mechanical 
resonator, (n), the photon count rate in each time bin then portrays the time-evolution of ( n)(t ) during the pulse 
on-state. All measurements presented herein were performed at a fridge base temperature of Tf = 11 mK. Further 
details about device fabrication and the experimental setup can be found in Refs. [13 EU and in the appendices. 

The signal-to-noise ratio (SNR) of this phonon counting method is determined by the sideband scattering rate 
70 m = 4«7q?t c /k (go is the optomechanical coupling rate, n c is the intracavity photon number, and k is the optical 
decay rate), the total system detection efficiency ( r ]), the dark count rate of the SPD (T dark ), and the residual 
transmission of the filters at the pump frequency relative to the peak transmission (A). A useful parameterization of 
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where n e is the optical decay rate into the detection channel. For the device under test we have k/2ty = 443 MHz, 
n e /2i r = 221 MHz, and go/2n = 710 kHz. A typical measured tjnep for our setup, taken using a comparable device at 
room temperature, is shown in Fig. as gray squares. Here, to obtain the lowest tinep optical pre-filtering is used to 
remove phase noise around the laser line and broadband spontaneous emission from the optical probe beam, details 
of which are given in App. [e] For sufficiently high probe power (n c > 40), u-nep falls below 10 —2 , enabling sensitive 
detection of the mechanical resonator deep in its quantum ground state. However, at sub-kelvin temperatures optical 
absorption heating produces a steady-state (n) > 1 for n c > 0.01 (red circles in Fig. 00 during continuous-wave (CW) 
optical excitation. In order to maintain the OMC in the mechanical ground state, the duty cycle of the pulse train 
must be kept sufficiently low, and the modulation depth sufficiently high, such that the mechanical mode thermalizes 
to the dilution refrigerator ambient bath between successive pulses. 

The CW behavior fits well to a thermal model consisting of the three processes illustrated in Fig [lj;: (i) the 
radiation pressure coupling at rate 70 M between the mechanical mode and the effective zero-temperature probe 
laser, (ii) coupling to an optical-absorption-induced hot phonon bath above the phononic shield bandgap, and (iii) 
coupling to the ambient fridge bath at an intrinsic rate 70 through the acoustic radiation shield. At the low intra¬ 
cavity photon numbers of these measurements, we believe the optical absorption heating is a result of excitation of 
electronic defect states at the silicon surfaces [23!, SI], and subsequent phonon-assisted relaxation of these states. As 
detailed in Ref. m, the resulting local hot-phonon bath occupancy (n p ) is found to scale as nJ A in steady-state, 
consistent with linear optical absorption and a cubic drop in the thermal conductance with temperature [25]. The 
corresponding coupling rate of the mechanical resonance to the high frequency hot-phonon bath ( 7 p ) is measured to 
scale as T p exp(—hui c /kBT p ) for low bath temperature (T p < 4 K), corresponding to inelastic phonon scattering with 
a quasi-equilibrium hot phonon bath above a cut-off phonon frequency of lo c /2tt « 35 GHz. As shown in Fig |TJ:1, 
extrapolation of this steady-state model for a fridge base temperature of Tf = 10 mK reveals a relevant pulse-off-state 
regime of n c < 10 -4 in which absorption heating effects are negligible. 

Figure [2]i shows the measured sideband photon count rate versus time with pulsed optical excitation, for both red- 
(A = w m ) and blue-detuned (A = — w m ) pumping. The on-state pulse amplitude in these measurements is n Ct0n = 45, 
corresponding to an optomechanical coupling rate of 7 om/27t = 205 kHz. For the achievable pulse extinction ratio of 
60 — 70 dB, this results in a residual off-state photon number of n c ,off < 4.5 x 10 -5 . Vertical dashed lines indicate 
the time bins corresponding to the start and stop of the pulse, determined from observing the rising and falling edges 
of the pulse when bypassing the cavity. The start and stop times are taken to be the time bins for which the pulse 
reaches 90% of its maximum value. Henceforth, the variable t refers to time relative to the synchronization signal 
generated by the pulse generator, while t pu ise refers to time relative to the start of the optical pulse occurring around 
t sa 1 fj, s. The pulse period in these measurements is fixed at T per = 5 ms. 

Throughout the pulse a pronounced asymmetry is observed between count rates for red-detuned versus blue-detuned 
pumping, which can be quantified by the asymmetry parameter £ = T_/r + — 1, where T± is the sideband photon 
count rate for a pump detuning A = ±oj m . While motional sideband asymmetry has previously been measured in a 
variety of opto- and electromechanical systems using linear detection schemes [20K221 [26] , the use of phonon counting 
techniques allows the observed asymmetry to be directly and unambiguously attributed to quantum fluctuations of 
the mechanical oscillator EH 123- This asymmetry, shown versus t pu ise in Fig. initially decreases with time before 
leveling off and beginning to increase for sufficiently long pulse times. The increase at later times can be ascribed to the 
effect of optomechanical backaction, which results in cooling or heating of the mechanical resonator for red- or blue- 
detuned pumps, respectively [23] ■ However, for t pu i se 7 qm ~ 750 ns the effects of backaction can be neglected and 
the phonon occupancy (n) may be assumed equal for both pump detunings. In this case, the asymmetry is simply 
related to the occupancy by £ = ( n ) _1 , and arises from the fundamental asymmetry between phonon absorption 
(r + a (n)) and emission (T_ a (n) + 1) processes [Ml [29]. Theoretical plots of the two contributions to £(f) are 
shown in Fig. §>■ The green dashed line shows the expected £(t) assuming that (n) (A = —w m ) = (n) (A = u m ) (no 
backaction), while the orange dash-dotted line shows £(f) in the absence of optical heating and in the case when 
rib ^ 1 such that £(0) =0 and asymmetry arises solely from backaction effects. From the asymmetry measured 
in the first 25 ns time bin, we extract a minimum phonon occupancy of (n) m ; n = 0.021 ± 0.007 (T m ; n ss 70 mK). 
Extrapolation to zero time yields a mode temperature closer to 40 mK, but still above the 11 mK of the mixing plate 
of the fridge, likely indicating a higher local chip temperature. 

Once (n) is determined in the initial time bin it may be used to convert measured count rates to phonon occupancies 
throughout the pulse for either pump detuning. The resulting calibrated occupancy versus t puIse is shown in Fig. [2};. 
A simple model of the heating dynamics would suggest an exponential increase towards the steady-state phonon 
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FIG. 2. a Total photon count rate versus time t for a red-detuned (red circles, A = w m ) and blue-detuned (blue squares, 
A = w m ). For both data sets n c ,on ~ 45 and Tper = 5 ms. Vertical dashed lines indicate the start and stop times of the pulses, 
b, Asymmetry £ versus time within the pulse t pu i se . The green dashed and orange dash-dotted lines show the theoretically 
expected £(t) in the absence of backaction and heating effects, respectively. Error bars show one standard deviation (s.d.) 
determined from the measured count rates, assuming Poissonian counting statistics, c, Calibrated phonon occupancy (n) 
versus t pu i se for A = u} m (red circles) and A = -w m (blue squares). The solid black lines show fits to a model including a slow 
exponential turn-on of the hot phonon bath. The inset shows detail for t pu i se < 750 ns on a linear scale for A = w m . Error 
bars show one s.d. determined from the measured count rates, assuming Poissonian counting statistics. For the blue-detuned 
data the error bars are smaller than the data markers. 


occupancy with total rate 7 = 7o + 7p + 7 om- However, as can be seen in the inset to Fig. ph, the curvature of the 
occupancy curve is positive for short times. This fact is inconsistent with the simple model described above, and 
is likely due to a finite equilibration time of the hot phonon bath. The data can be fit well for both detunings by 
using a simple model which assumes that a small fraction of the hot phonon bath turns on effectively instantaneously 
(i.e. very fast relative to the length of a single time bin) while the remainder has a slow exponential increase to its 
steady-state value. Thus, the effective rate equation for the phonon occupancy is given by 

(n) = - 7 (n) + 7 p ™p (l - <5be -TstpulM ), (2) 

where <5b and 7s are the slow growing fraction of the hot phonon bath and the corresponding turn-on rate, respectively, 
and the steady-state hot phonon bath parameters 7 P and n p can be determined by fitting the CW (n) versus n c data 
shown in Fig. [lji. Fitting the pulsed occupancy curve yields <5b = 0.79 ± 0.08 and 7s/27r = 215 ± 29 kHz, indicating 
that the heating occurs slowly enough as to be manageable during coherent quantum operations as determined in the 
remaining discussion. 

During the off-state of the pulse, optical heating of the mechanical resonator should be negligible and the phonon 
occupancy should cool at the intrinsic damping rate 70 . Thus, the initial and final occupancies during the pulse ((n)i 
and (n)f, respectively) should obey the relation (n )i = e _7oTper (n)f, assuming the pulse period T per is much larger 















5 



FIG. 3. Ratio of initial/final phonon number ( n)i/(n)t versus pulse period Tper- The solid red curve shows a fit to exponential 
decay. Error bars show one s.d. determined from the measured count rates, assuming Poissonian counting statistics The inset 
shows (n) versus t pu iae f° r the different values of T pel . 
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FIG. 4. a, Maximum phonon occupancy during a pulse, (n) ma x, versus pulse period Tper and pulse width T pu i se . The white 
contour delineates the region where the effective cooperativity C e a > 1 throughout the pulse, b, C e ff versus time within the 
pulse f pulse for Tper = 5 ms, Tpuise = 3 fis. 


than the pulse width T pu i se . The ratio (n)i/(n)f, shown in Fig. [3^, versus Tper, displays the expected exponential decay 
with 7o/2-7t = 328 ± 14 Hz and a corresponding intrinsic mechanical quality factor of Q m ~ 1-7 x 10'. This decay rate 
corresponds well with the value inferred from previous CW measurements of occupancy at mK temperatures |18] . From 
our measurement of 70 and (n) m i n we estimate a thermal decoherence time of r t h = ( 70(1 + (ti)min )) -1 = 475 ± 21 /is. 

While the low thermal occupancy and long coherence time measured here are promising, the utility of cavity 
optomechanical systems for performing coherent quantum operations between the optical and mechanical degrees 
of freedom is ultimately predicated upon the ability to simultaneously achieve (n) -C 1 and large cooperativity 
C = 70 M/ 7 bi where 7 b = 70 + 7 P is the total coupling rate between the mechanical resonator and its thermal bath. In 
the specific example of optomechanically mediated coherent transfer of photons between optical and superconducting 
microwave resonators M, the relevant figure of merit is the effective cooperativity C e s = C/n\>, where rib is the 
effective bath occupancy defined such that 7 brib = 7 o^o + 7 P ^ P i which must be much larger than unity in order to 
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achieve low-noise photon conversion at the single quantum level uni Eg. Using the measured values of 70 and 73 , 
we can calculate the maximum phonon occupancy achieved during the pulse, (n) max , for a given T pu i se and T per and 
A = w m (Fig. [4^,), as well as C eff as a function of t pu i se (Fig. |4jo). 

Though optical heating of our devices prevents us from reaching the C e g > 1 regime using a CW pump, due to 
the slow turn-on of the hot phonon bath observed in this work we find that (n) < 1 and C e s > 1 can be maintained 
throughout the entire pulse period for T pu i se < 300 ns at pulse rates approaching 1 MHz. Furthermore, we find 
C e ff 1 during the initial 100 ns of the pulse and reaches values as large as C e ff ~ 40. An analysis of the effects 
of optical heating on pulsed phonon addition and subtraction |5J for the creation of non-Gaussian mechanical states 
are presented in App. E wherein we calculate the fidelity of generating a mechanical Fock state |7] to be F > 0.8 
for pulse rates up to 1 kHz. These calculations demonstrate that optomechanical applications in the quantum regime 
are feasible with our current structure for reasonable pulse parameters without further mitigation of optical heating 
effects. 
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FIG. 5. Experimental setup. Full setup for performing pulsed phonon counting. A-meter: wavemeter, DeMux: C-band 
optical demultiplexer, (j>- m: electro-optic phase modulator, SW: optical switch, a-m: high-extinction electro-optic amplitude 
modulator, switch: high-speed electro-optic switch, VOA: variable optical attenuator, EDFA: erbium-doped fiber amplifier, 
PD: high-speed photodetector, RSA: real-time spectrum analyzer, VNA: vector network analyzer, SPD: single photon detector, 
TCSPC: time-correlated single photon counting, FPC: fiber polarization controller.. 
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Appendix A: Device Fabrication 


The devices are fabricated from a silicon-on-insulator (SOI) wafer (SOITEC, 220 nm device layer, 3 /mi buried oxide) 
using electron beam lithography followed by reactive ion etching (RIE/ICP). The Si device layer is then masked using 
ProTEK PSB photoresist to define a mesa region of the chip to which a tapered lensed fiber can access. Outside of the 
protected mesa region, the buried oxide is removed with a plasma etch and a trench is formed in the underlying silicon 
substrate using tetramethylammonium hydroxide (TMAH). The devices are then released in hydrofluoric acid (49 % 
aqueous HF solution) and cleaned in a piranha solution (3-to-l H 2 S 04 :H 202 ) before a final hydrogen termination in 
diluted HF. In fabrication, arrays of the nominal design are scaled by ±2 % to account for frequency shifts due to 
fabrication imperfections and disorder. 












































































Appendix B: Experimental Setup 


The full measurement setup is shown in Fig. [5] A fiber-coupled, wavelength tunable external cavity diode laser is 
used as the light source, with a small portion (~ 1%) of the laser output sent to a wavemeter (A-nreter) for frequency 
stabilization. The remaining laser power is sent through both a C-band optical demultiplexer (DeMux) to reject 
broadband spontaneous emission from the laser, and a high-finesse tunable filter to remove laser phase noise at the 
mechanical frequency, which can both contribute to excess pump transmission noise on the single photon detector [Hi- 
After pre-filtering, the laser pump is sent through an electro-optic phase modulator (</>-m) which may be driven by 
either an RF signal generator, to generate optical sidebands for locking the filter cavities, or a vector network analyzer 
(VNA) for obtaining the amplitude and phase response of the optical cavity. 

To generate pulses, an optical switch (SW1) allows the pump tone to be sent through a series of electro-optic 
amplitude modulators for the generation of optical pulses. The first is a fast high-extinction modulator (a-m) with a 
rise and fall time of ~ 25 ns. Though in principle the high-extinction modulator can provide 50 — 60 dB extinction on 
its own, its transmission level is much less stable and difficult to lock at the maximum extinction point. Thus, we lock 
the modulator at ~ 30 dB of extinction and use two electro-optic switches, each providing ~ 18 — 20 dB extinction, 
to fully extinguish the pump. While these switches are much slower (~ 200 ns rise time, ~ 30 /is fall time), the extra 
switching time should have a negligible impact as it is much smaller than the pulse periods used here (~ 1 ms or 
more). The modulators are collectively driven by a variable delay electrical pulse generator which also provides a 
synchronization pulse to the single photon counting electronics. 

A variable optical attenuator (VOA) controls the power input to the cavity, after which an optical circulator routes 
the laser light to a lensed fiber tip inside the dilution refrigerator for end-fire coupling to the device. Subsequently, 
the cavity reflection can be switched (SW3) into one of two detection setups. The first setup sends the signal 
through an erbium-doped fiber amplifier (EDFA) followed by a high-speed photodetector (PD). The resulting amplified 
photocurrent may be directed to a real-time spectrum analyzer (RSA) in order to measure the optical noise power 
spectral density (NPSD) for mechanical characterization or to the VNA which is used in conjunction with the phase 
modulator to measure the full complex response of the optical cavity for purposes of optical characterization as 
described below. The second detection path sends the reflected signal through three additional tunable filters in order 
to reject the pump frequency, and then back into the dilution refrigerator where it is detected by a superconducting 
single photon detector (SPD). The output of this detector is sent to a time-correlated single photon counting (TCSPC) 
module to build up a histogram of photon count events as a function of time relative to the sync pulse received from 
the pulse generator. 

The SPDs used in this work are amorphous WSi-based superconducting nanowire single-photon detectors developed 
in collaboration between the Jet Propulsion Laboratory and NIST. The SPDs are design for high-efficiency detection 
of individual photons in the wavelength range A = 1520 — 1610 nm with maximum count rates of about 2.5 x 10 7 counts 
per second (c.p.s) EDI- The SPDs are mounted on the still stage of the dilution refrigerator at ~ 700 mK. Single-mode 
optical fibers (Corning SMF-28) are passed into the refrigerator through vacuum feedthroughs and coupled to the 
SPDs via a fiber sleeve attached to each SPD mount. Proper alignment of the incoming fiber with the 15 /mi x 15 /tm 
square area of the SPD nanowire is ensured by a self-aligned mounting system incorporated into the design of the 
SPD [3H] . The radio-frequency output of each SPD is amplified by a cold-amplifier mounted on the 50 K stage of 
the refrigerator as well as a room-temperature amplifier, then read out by a triggered PicoQuant PicoHarp 300 time- 
correlated single photon counting module. By systematically isolating the input optical fiber from environmental light 
sources and filtering out long wavelength blackbody radiation inside the fridge we have achieved dark count rates of 
~ 4 (c.p.s.). At just below the switching current of the detectors, we have measured a peak detection efficiency of 
Vspb = 68%, with < 20% variability depending on photon polarization. 

The tunable filters used for both pre-filtering the pump and filtering the cavity reflection are commercially available, 
piezo-tunable Fabry-Perot filters (Micron Optics, FFP-TF2), all with a ~ 50 MHz bandwidth, a free-spectral range 
of ~ 20 GHz, and a tuning voltage of < 18 V per free-spectral range. The filters each offer roughly 40 dB of pump 
suppression (relative to the peak transmission) measured at a filter-pump varies by 1 — 2 dB from filter to filter. When 
locking the post-cavity filters a switch is used to bypass the cavity (SW2), as a relatively large amount of CW power 
is used during the locking procedure and we would like to avoid sending large amounts of power into the cavity unless 
necessary. More importantly, the SPD is also bypassed (SW4), as allowing too much power to reach the SPD will 
saturate it resulting in significantly elevated dark counts (Fdark ~ 500 — 1000 c.p.s.) for 1 — 2 minutes after the signal 
is turned off. Once the switches are set an RF signal generator is used to drive the phase modulator at w m , producing 
an optical sideband which is aligned with the cavity (and motional sideband) frequency. To lock the filter chain, a 
sinusoidal voltage (with an initial range of ±10 V) is used to dither each filter while monitoring its transmission. 
The offsets of the sinusoids are then adjusted, and their amplitudes reduced, to maximize transmission of the desired 
sideband, fixing the voltages once all three filters are well aligned. Over time the filters will drift and the locking 
procedure will need to be repeated, though in subsequent re-locking attempts a much smaller dithering amplitude 
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of ~ 1 V is sufficient. As the piezo elements controlling each filter have a finite settling time, the filters will drift 
much more rapidly following the initial locking attempt (during which large voltage shifts are applied) than following 
subsequent re-locks. After several re-locks, the filters will typically become stable enough that the total transmission 
at the sideband frequency changes by < 5% over several minutes. At this point the phase modulator is turned off, the 
pulse generator, cavity and SPD are switched back into the optical train and we begin accumulating pulsed data for 
2 minutes before re-locking the filters. The total filter transmission is recorded at the end of a locking procedure and 
again prior to re-locking, and if the transmission has shifted by more than a few percent the previous pulsed dataset 
is discarded. All the time-resolved data shown in this chapter are taken in this manner, with all 2 minutes datasets 
averaged together to produce the final pulsed histogram. 


Appendix C: Device Characterization 

The measurements presented in this work rely on an accurate calibration of the optomechanical damping rate 7om ; 
which depends on the vacuum optomechanical coupling rate go, the total optical decay rate k, and the intracavity 
photon number n c . The photon number for a given power and detuning in turn depends on the single pass hber- 
to-waveguide coupling efficiency g cp i and waveguide-cavity coupling efficiency rj K = n e /n (re e is the waveguide-cavity 
coupling rate). The fiber collection efficiency is determined by measuring the calibrated reflection level far-off resonance 
with the optical cavity on the optical power meter, and is found to be r] C p\ = 0.68 for the device measured here (total 
device reflection efficiency of ~ 0.46%). To measure k and r] K , the laser is placed off-resonance from the cavity and 
the VNA is used to drive the phase modulator and sweep an optical sideband across the cavity. By detecting the 
reflection on a high-speed photodiode connected to the VNA input we obtain the amplitude and phase response of 
the cavity from which can extract k/2tt = 443 MHz (corresponding to an optical quality factor of Q a « 4.4 x 10 5 ), 
and r] K = 0.5. With these three parameters measured, it is possible to determine n c for an arbitrary input power to 
the cavity. 

To characterize the acoustic resonance, the EDFA is used to amplify the cavity reflection so that the optical 
noise floor overcomes the photodetector’s electronic noise and the optical NPSD is measured on the RSA, where a 
Lorentzian response due to transduction of the acoustic thermal Brownian motion can be observed at the acoustic 
resonance frequency w m /27r = 5.6 GHz. For a red- or blue detuned pump laser (A = to c — ui\ = ±w m ) the linewidth 
of this transduced noise peak is 7 = 7 b ± 7om, where 7 b is the bare coupling rate of the mechanical resonator to 
its bath and 70 M = ^ 9 qTi c /k. Thus, by averaging the observed linewidth for red- and blue-detuned pumps we can 
determine 7b and thus extract the excess optomechanically induced damping rate 70 M as a function of n c . With n c 
and k determined, a linear fit of 70 M reveals a vacuum optomechanical coupling rate of <?o/27t = 710 kHz. 


Appendix D: Phonon Counting 


In the linearized sideband resolved regime with a red-detuned pump (A = w m ) the total reflected cavity amplitude 
will be sum of the coherent pump reflection a out and a fluctuation term given, in the Fourier domain and in a frame 
rotating at the pump frequency, by m 

a 0 utMU=a; m ~ r(w)a in (w) + n(w)ai(w) + s(uj)b in (uj ), (Dl) 


where di„ and ch are the pump noise and intrinsic cavity noise, respectively (assumed for the moment to be vacuum 
noise), b ln is the noise operator for the mechanical bath, which has the correlation function (i>l n (u})b^ (co 1 )^ = no8[u>-\- 
io') (no is the bath occupancy), and the effective scattering matrix elements are given by 


r(w) = 
n(u;±) 


2ft e _j_ 70M^e 1 

K K ±i (w m =F w) + 7/2 ’ 

j_ V K i K e ( _ 7QM _ 2 "\ 

K \±i (w m I«)+ 7/2 ) ’ 


s(w) 


V K ±* (Um T w) + 7/2 


(D2) 

(D3) 

(D4) 


As the full reflected signal from the cavity includes the reflected pump amplitude, which is many orders of magnitude 
larger than the sideband signal, we must first filter the cavity output to reject the pump. Provided the filter is 
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sufficiently high-finesse, it can be modeled in the frequency domain by a single Lorentzian function 

«f/2 


F{w; Wf) = t 


i(ui{ — to) + k { /2 ’ 


(D5) 


where k,{ and uif are the resonant frequency of the filter, respectively. Now, explicitly considering the case of red- 
detuned driving in the sideband-resolved regime, the filtered cavity output will be the product of with the 

frequency domain output of the cavity. As the resonantly enhanced anti-Stokes sideband photons will be detuned by 
w m from the pump, we choose ojf = w m . Using Eq. |D1[ we obtain 

dfiit(w) = F(w;w m ) ^a out (5(w)+r(w;+)d in (w)-t-n(w;+)ai(a;)+s(w;+) 6 in (w)^ . (D 6 ) 

Performing photon counting on the filtered output then results in an average count rate of 

r(f) = 

oo 

duj I du' e l ^ +u>,)t 

J —OO ' 

1 / K r°° 

7 T ( l^(°; w m)| 2 |aout| 2 + —■tom / dw|F(w;w m )| 2 S hh ( uj -, ( n )) 

2tt V K J ~oo 


mt v 

Y /“OO /“OO 

2n 


27r 

“4|a 2 ut | + ^ 70 M(n), 


(D7) 


where A = A-|F(0;w m )| 2 is the pump attenuation factor, and S’ bb (w; (n)) is the phonon spectral density 


5 bb (w) = 


q(n) 


(w m + w) 2 + (7/2) 2 ’ 


(D8) 


Note that we have assumed a filter bandwidth Kf 7 , allowing us to approximate |F(w;w m )| 2 ~ |F(w m ;o ; m )| 2 = 1 
inside the integral over Sbb- A similar analysis for blue detuning (where Wf = — w m to filter the Stokes sideband) 
yields a comparable result, with (n) —> (n) +1. Note that this formula assumes that Si n and d; are vacuum noise, and 
thus do not contribute to counting of real photons. 

The total count rate, including noise of the photon counter and reduction of the pump and signal due to measurement 
inefficiency, is given for red-detuning by 


Ttot — I 'dark A IV 


+ V— 70M(n), 
Hi 


(D9) 


where Tdark is the dark count rate of the photon detector, r pump = r/A\a ou t\ 2 and rj is the total measurement 
efficiency. These expressions can be used to perform thermometry in a similar fashion to linear detection, either by 
calibrating the cavity parameters and total measurement efficiency or by measuring the asymmetry between the red- 
and blue-detuned count rates. 

To assess the sensitivity of this counting scheme, it is convenient to express the measurement noise floor Tdark+rpump 
in terms of an equivalent number of mechanical quanta (that is the mechanical occupancy (n) that would be needed 
to yield a signal-to-noise of one). This noise-equivalent phonon number is obtained by dividing the total noise floor 
by the per-phonon count rate Tsb,o = v( K e /yielding 


«nep = 


r dark + T 


pump 


Tsb.o 


(DIO) 


For a highly sideband-resolved system, the reflected pump in the case of A = ±w m will be approximately given by 
owt ~ ol j n . This in turn can be expressed in terms of the intracavity photon number as |Ck ou t1 2 ~ w^n c /K e . Thus, 
^nep as a function of n c is given by 


k 2 T 


dark 


A 



2 


%Ep(n c ) 


4 rjKe.gfin, 


2 ^e ^0 


(Dll) 
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Appendix E: Effects of Technical Laser Noise 

The mechanical frequency of the nanobeam used in this experiment (w m / 2ir = 5.6 GHz) raises concerns about the 
effects of laser phase noise on the measurements, as the laser used in this experiment has previously been observed 
to possess a prominent phase noise peak at 5 GHz m- In addition to phase noise, most diode lasers typically have 
a small amount of broadband spontaneous emission. While this additional noise is orders of magnitude weaker than 
the laser tone itself, it exists outside the wavelength region (A ~ 1520 — 1570 nm) where the filters are guaranteed to 
be high finesse, and thus can be transmitted with high efficiency to the SPDs. 

Phase noise in particular is worrisome as it can not only lead to an excess noise floor but also to real heating of the 
mechanics and systematic errors in thermometry due to noise squashing/anti-squashing |321133] . Phase noise can be 
accounted for by assuming a total pump noise operator given by [3D 

<2in,tot(0 = Uin(t) T (El) 

where &i n (f) still represents vacuum noise and 

Hit) = i\a in \</){t), (E2) 

where |a; n | is the amplitude of the pump and is the stochastic phase of the pump, assumed to be a real, stationary 
Gaussian process with zero mean. The phase noise is assumed to be delta-correlated in the frequency domain, such 
that 


(</>(w)</>(u/)) = S^{u)8{u + u/), (E3) 

where the expectation value here corresponds to an ensemble average. The phase noise input to the cavity then has 
the correlation function 


(|aj,(w)a0(w ')^ = S aa (uj)5(ijj + u'), (E4) 

where S aa (u>) = \a\ 2 S^(u). In a sideband-resolved system, for either red or blue detuning A = ±w m , we find that 
the presence of phase noise heats the mechanical resonator, with an additional added phonon occupancy given by [5T1 

/ \ I 70M /rir\ 

\n)4,\A=±u m = -(E5) 

k 7 

where we have defined = S aa (u> m ), and where we have assumed that S aa (u> ) is sufficiently slow-varying in the 
vicinity of ui = (specifically for \ui — w m | < 7 ) that we may approximate S aa (oj) = S aa (io m ). 

Including this additional noise term in the analysis of the previous section yields a phase noise contribution to the 
total photon count rate of 


/ OO 7 

—S aa (uj)\F(uj, ±w m )| 2 |r(w; ±)| 2 

-OO 


If we assume that S aa (uj) is slowly-varying in frequency for | u> — w m | < 7 and that Kf 7 , this simplifies to 

2 


r 


01 A. 




=±^ m =W4> j 


1 - 


2n e 


«e7QM / 7QM^e _j_ A _ 2fC P 
K \ 7 K 


(E 6 ) 


(E7) 


Using the fact that = \a\ 2 S^iuiro) ~ w^ 1 n c 5 , 00 (a; m )/ K e , we obtain the contribution of phase noise to the noise- 
equivalent phonon number 


«NEP,0|A=±w m — 


2K e go 




Kf 


1 - 


2n f 


K e70M / 70M^e A _ 2n e 
K \ 7 K 


(E 8 ) 


Like the pump-bleed through, phase noise leads to a constant contribution to Unep, and leads to squashing or anti¬ 
squashing of the noise depending on detuning and cooperativity, similar to heterodyne detection. Note, however, that 
in the case K e /n = 0.5 the contribution of phase noise will not depend upon detuning. Thus, even in the presence of 
large phase noise it is possible to avoid detuning dependent noise squashing/anti-squashing, though one will still have 
a large overall phase noise floor. 
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It is useful for characterization purposes to calculate the phase-noise contribution to the observed count rates when 
the laser is far-detuned from the cavity resonance (| A| w m ). Assuming that the laser-filter detuning is kept constant 
at ±w m , the phase-noise count rate in this case is just 


/ oo 

— Sj M (w)|.F(u;,±u; m )|" 


K[ 


with a corresponding noise-equivalent phonon number 


^NEP,0||A|>tj m 


/ UJ m K 

\4 n e go 


2 

^{S(j)(j) (uj m) • 


(E9) 


(E10) 


To get rid of the excess noise, we insert both the bandpass filter (for filtering spontaneous emission) and a tunable 
high-finesse filter (for filtering phase noise) immediately after the pump laser output as shown in Fig. [5j enabling us 
to reach «nep 1 using a three-filter phonon counting setup as shown in Fig. Id of the main text. A conservative 
estimate of the residual phase noise can be made by assuming that the limi ting value of unep ~4x 10~ 3 is entirely 
due to phase noise (i.e. perfect filtering of the pump tone). Using Eq. E10 we find S^(w m ) <8x 10 -19 Hz -1 . For 
the pump power during the on-state of the pulse (n Cj0n ~ 45), the corresponding effective phase noise occupancy is 
« 3.2 x 10~ 5 , which has a negligible effect on the measurements in this work. 


Appendix F: Heating model 


The simplest thermal model assumes that the optically induced bath turns on instantaneously when the pulse is 
in the on-state. The corresponding rate equation for the phonon occupancy (n), for red- (A = w m ) and blue-detuned 
(A = — w m ) pulses during the on-state is thus 


(n) 


-7(n) + 7 P n p + 7 0 n 0 + - (1 T 1) 70M, 


(FI) 


where 7 = 7 o + 7 P ± 7om, 7 p and n p are the coupling rate and occupancy of the hot phonon bath, 70 and no are the 
coupling rate and occupancy of the ambient fridge bath, and the extra factor of 70 M for a blue-detuned pump accounts 
for the possibility of spontaneous emission into the mechanical resonator due to the optomechanical interaction. This 
rate equation has the simple solution 


(n)(t) = (n)(t 0 )e 7 ‘+ n H (1-e 7 ‘) , (F2) 

«h = 7 _1 [bp^p + 7o«o + * (1 =F 1) 70M j ■ (F3) 


where to is the start time of the pulse, and to < t < to + T pu i se . 

In principle, 70 M can be determined independently as described above, 70 can be determined by fitting the occu¬ 
pancy decay during the pulse off-state (Fig. 3 in the main text), while 7 p and n p can be subsequently determined 
by fitting the steady-state occupancy curve shown in Fig. lc of the main text. However, using these independently 
determined values in a fit to the red- and blue-detuned data shown in Fig. 2a of the main text results in a poor fit 
and inconsistent results. In particular, the apparent heating rate 7 is much smaller than expected for a red-detuned 
pulse and larger than expected for a blue-detuned pulse. This, along with the “kink” in the heating curve shown in 
the inset of Fig. 2c of the main text, suggests additional complexity in the heating dynamics. 

The simplest addition to the heating model is to assume a finite time for the hot phonon bath to come into 
equilibrium, which is approximated by allowing a fraction of the hot phonon bath occupancy to turn on exponentially 
over time. Thus, the phenomenological rate equation becomes 

(n) = -7 (n) + 7 p n p (l - 6 b e _7st ) + 7 0 n 0 + ^ (1 T 1) 70 M, (F4) 

where 5 b is the slow growing fraction of n p and 73 the turn-on rate. Strictly speaking 7 p should depend on the 
phonon distribution of the hot phonon bath, and thus would be expected to be time-dependent in this model as well. 
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However, the resulting rate equation becomes intractable in such a case and the effects should be negligible in the 
regime 7 om ^ 7 P , so we approximate 7 P equal to its steady-state value. This modified rate equation has the solution 

(' n){t) = (n)(f 0 )e _7t + n H (l - e -7 *) + n s (e _7st - e _7 ‘) , n s = 7pn P^ b , (F5) 

7s 7 


which is used to obtain the fit shown in Fig. 2c of the main text with (n)(to), 7 s and <5b as free parameters. 

During the off-state of the pulse (to + Tpuise < t < to + T per ), the resonator will simply cool towards the ambient 
fridge occupancy no at the intrinsic damping rate 70 . Using the initial condition that (n)(to) = no for the first pulse 
(to = 0 ), and iterating many pulses we find that in the steady-state the initial phonon occupancy during a pulse 
(assuming T per > T pulse ) is 


\,. , _ n 0 (1 ■ 


(F 6 ) 


Thus, once n 0 , 7 s and (5b are determined by fitting the occupancy curves, we may use Eqs. |F5| and |F6| to determine 
the occupancy throughout the pulse for arbitrary T per and T pu i se . This allows us to determine the maximum attained 
phonon occupancy as a function of pulse parameters, as shown in Fig. 4a of the main text. 


Appendix G: Phonon addition/subtraction fidelity 

In general, the optomechanical interaction allows for the creation of a variety of non-classical mechanical states 
via phonon subtraction or addition using appropriate red- or blue-detuned pulses [5]- In particular, recent proposals 
have explicitly shown how to herald single phonon Fock states [7j and entangled mechanical states [ 6 ] via application 
of a short blue-detuned pulse to an optomechanical system in its motional quantum ground state and subsequent 
detection of the emitted sideband photon. However, the fidelity of the generated states in our system is in principle 
degraded by the presence of optical heating during the pulse, which we analyze to lowest order here. 

In a frame rotating at the pump frequency, the full Heisenberg-Langevin equations for the optomechanical system 
are |5L 


a = — ^iA + a + ig 0 (b + 6^ a + y^din, 
b = - (iu m + y) b + ig 0 a^a + y/jbb in , 


(Gl) 

(G2) 


where a, b are the photon and phonon annihilation operators, respectively, and di n , b m are quantum noise operators 
corresponding to the optical and mechanical baths. 

We linearize about a classical steady-state by displacing a —> a + a, where \a\ 2 = n c . For concreteness we will 
consider the case A = w c — ui\ = — oj m (blue-detuned pump). Moving into a frame rotating at the mechanical frequency 
(i.e. a —> ae IWmt , and so on for all operators), and making the rotating wave approximation, valid in the weak coupling 
(Soy^c k) and sideband-resolved (ft/w m “C 1 ) limit, we obtain 


a = ^ 2 ® + V 

h=-^b + iGal + ^b in , 


(G3) 

(G4) 


where G = anc ^ ®in> b\ n are the usual noise operators multiplied by e The noise operators obey the 

following commutation and correlation relations 


a in {t),al n (t') = b in {t),bl n (t') =8{t-t’) 

= S(t — t') 

^in (t)bin{t')) = n h (t)6{t - t'), (bin{t)b[ n (t')) = ( n b (t ) + 1) 5(t - t'). 


(G5) 

(G 6 ) 

(G7) 


Since we are working in the weak-coupling limit (G <C n) we may use the adiabatic solution for a (i.e. a « 0). 
Moreover, we wish to include the effects of mechanical noise to lowest order. Considering the effect of a short optical 
pulse of duration r, we consider the case 7b t, 7b / Q r ds rib(s) <C 1 , as well as 7 omt 1 and 70M 5§> 7b, where here 
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7om = 4 G 2 /k refers to the magnitude of the optomechanical damping rate (the sign will be explicitly incorporated 
for simplicity, since we’re only considering blue-detuning). Under these assumptions, as in Ref. [34], the mechanical 
noise term will be retained in the mechanical equation of motion only, since in the optical equations of motion it 
will acquire an extra factor of 7om t and thus can be neglected to lowest order. Furthermore, the (7b/2) 6 term will 
be neglected in favor of the much larger (7om/2) 6 term. With these approximations, we arrive at the approximate 
equations 


TOM 


6 f + -din, 


7om 


b + ^VTOMdjL + VTbdin. 


(G 8 ) 

(G9) 


Formally integrating the equation of motion for b over the duration of the pulse while using the input-output relation 
for the optical cavity d ou t = — di n + ^/Kd, we arrive at the approximate equations 


t t 

^OUt(^) = ttin (t)+i«e^St( 0 ) + 7 oMe^ ds 

Jo 


b(t) = e 2 6 ( 0 ) + iyj 7 QMe 2 


ds 


„f / \ .- 'OM 

hn( s ) + \/ 7 b e 2 


— TQM s „ , . 

e 2 a in (s) 


--TQM» y , . 

ds e 3 Oin(s). 


We now introduce the following temporal modes [34. 

Bin = 6 ( 0 ), -Bout = 6 (r), 

-din = 


7om 


-dout — 


1 — g—T om’" 
7QM 

gTOMT _ ^ 
'TOM 7 ' / 


-TQM- „ , . 

ds e 2 a; n (s), 


TQM- „ 

ds e 2 a out (s) 


ds > 


bin(s ), 


/ 0 


F = 

which allows us to rewrite the input/output equations at the end of the pulse it = r) 

iout = e 2 ^i in + iVe JOMT - 1 Bj n , 

B out = An + fVeTOMr _ lit + p_ 


(G10) 

(Gil) 

(G12) 

(G13) 

(G14) 

(G15) 


(G16) 

(G17) 


Note that [An, An] = [Amt, ALt] = [An, An] = 1 as expected. As for the noise operator F, the commutator and 
correlation functions are defined, using the known properties of to lowest order as 

[F, F f ] w 7 b r + 0(7oM7bT 2 ), (G18) 

<F f f} = 7 b e 7 ° MT J ds e- 7 ° MS ?r b (s), (fF^ = 7b r + <F* A} • (G19) 

We now seek a unitary propagator U such that A out = WA- ln U and B ou t = WB ln U. Note that if we define 
e 7 oM-r /2 _ C osh(r), •v/e 7 ° MT — 1 = sinh(r), then the input/output equations have the form of two-mode squeezing 
between i; n and B; n in the absence of mechanical noise F. Thus, we want U = B m A> where 


B 0 = e -(A n B i „+Tt n st 1 ) ! 

(G20) 

AnAnAn = An + A 

(G21) 

[An, An] = 0. 

(G22) 


We find that a beam-splitter type interaction between B; n and F satisfies this condition to lowest order. That is, 
if An = . then using the above commutation relations for F we can show, 

ALAnAn = COS (x/ 7 b r) An + (7bT ) _1/2 sin (^yjyr) F 
~ An + F. 


(G23) 





























15 


Given this propagator U which yields the approximate system evolution over the pulse interval r, and assuming 
an initial density matrix p o = |0a) (Oa 8 Pb,o 8 Pf, where A, B, and F refer to the optical, mechanical, and bath 
subsystems respectively, the density matrix for the mechanical resonator, p b, conditioned upon detection of a single 
photon, is given by 


TrA.F (|1 a> <1a| <8 Ib 8 IfU PoA) 
Tr (| 1a) (1a | 8 Ib 8 IfPo) 


= F " 1 Tr F (An ( 1 A | U 0 P 0 U& |U) AL) , 

(G24) 

where P is the overall probability of photodetection, given by 


F = Tr B ((1 A | A |0 A ) (0 a | 8 Pb.oAJ |U)) ■ 

(G25) 

Denoting ps, ideal = P 1 AnPB,oAn as the ideal conditional matrix in the absence of mechanical dissipation, and 
expanding U to lowest order in 7 omG 7 bT and ^F^F 1 ), we find 

F « 70MT (AnAn) , 

(G26) 


and 


PB ~ (l — (Af)) PB, ideal + (.Af) AnPB.idealAn + ^7b T + (Af)) AnPB,ideal An 

— 2 (^ bT + ^(Af)) ^AnAnPB,ideal + PB.idealAnAn) • 


(G27) 


To evaluate the fidelity of the phonon addition subtraction, we use the definition of fidelity between two quantum 
states pi and pa, given by [35] F = Tr (\Jp^ 2 piP \ 2 \, which is straightforward to evaluate for an initial thermal 


mechanical state, as both pi and p 2 will be diagonal in the Fock state basis. 

Note that the only substantive difference in the case of phonon subtraction (Delta = w m , red-detuned pumping) is 
that the propagator U 0 will have the form of a beam-splitter interaction rather than two mode squeezing. This will 
only change P and ps, ideal, which are now given by 

(G28) 
(G29) 


P ~ 70MT ^/ri^in) > 
PB, ideal ^inPB,0-^in’ 


while the definition of pb in terms of ideal is unchanged assuming that the bath heating is only dependent on pump 
power and not pump detuning. 

To evaluate the fidelity of the creation of a single phonon Fock state [7] in our system as a function of pulse width 
and period for blue-detuned puls es, w e assume the initial state of the mechanics to be a thermal state with average 
phonon number (n) given by Eq. 


F 6 


while (Af) can be evaluated using the fit parameters of the effective time- 


dependent hot phonon bath. We may then calculate the conditional density matrix via Eq. |G27| and the fidelity by 
setting pi = |Ib) (Ib| and p 2 = Pb- Considering a range of pulse widths 0 < T pu i se < 100 ns, which is sufficient to 
guarantee 7 oMFpuise 1 , we find that for T per > 1 ms there is no appreciable loss in fidelity due to long term heating 
causing an increase in ( n)(to ), and the primary loss in fidelity is due to transient heating during the pulse. Thus, 
we set the F per = 1 ms and calculate the fidelity of Fock state generation as a function of pulse width, displayed in 
Fig#- For short (F pu i se < 10 ns) pulse widths, the fidelity approaches 98.5%, and remains above 80% for pulse widths 
approaching 100 ns. However, it is equally important to quantify the expected time to herald such a state (i.e. the 
expected time before detection of a sideband photon), which in this case is given by Fpock = Fper/(p70MFpuise(n)(fo)), 
where p is the total detection efficiency of the phonon counting measurement. This is shown versus F pu i se in Fig. 
for two cases. The solid line is calculated the actual measured detection efficiency of p ss 0.3%, while the dashed 
line is calculated using a realistic estimate of the ideal efficiency of our current measurement setup, p ss 5.5%. This 
latter figure is calculated using the measured efficiencies of the SPD (pspd = 68 %), the fiber-to-waveguide coupling 
efficiency (p cp i = 68 %), and the waveguide-cavity coupling efficiency (p K = 50%), and adding an additional 2 dB 
insertion loss per filter and 0.5 dB for the insertion loss of the optical circulator, which correspond to the highest 
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FIG. 6. Fock state generation, a, Fidelity F of the generation of a single phonon Fock state versus pulse width T pu i se 
(pulse period T pel = 1 ms), b, Average time required to herald a Fock state TFock versus pulse width T pu i se . The solid line is 
calculated using the total detection efficiency measured in this work, while the dashed line uses the estimated ideal detection 
efficiency. 


efficiencies measured for these components in our lab. As shown in Fig. [6j even in the idealized case the expected time 
for the generation of a Fock state is > 100 ms, which is much longer than the lifetime of the mechanical state. Thus, 
while it is feasible to use our current systems for the heralded, high-fidelity generation of non-classical mechanical 
states, it it still necessary to reduce heating effects (and thus allow for shorter pulse periods) in order to utilize this 
procedure for useful quantum information processing tasks (e.g. scalable entanglement distribution via the DLCZ 
protocol EM). 





